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We discuss the possibilities of detecting radio pulses 
from high energy showers in ice, such as those produced 
by PeV and EeV neutrino interactions. It is shown that 
the rich radiation pattern structure in the 100 MHz to 
few GHz allows the separation of electromagnetic show- 
ers induced by photons or electrons above 100 PeV from 
those induced by hadrons. This opens up the possibil- 
ity of measuring the energy fraction transmitted to the 
electron in a charged current electron neutrino interac- 
tion with adequate sampling of the angular distribution 
of the signal. The radio technique has the potential to 
complement conventional high energy neutrino detectors 
with flavor information. 
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High energy neutrino detection is one of the experi- 
mental challenges for the next decade with efforts under 
way to construct large Cherenkov detectors arrays under 
water or ice Q. The size of these detectors must be in 
the scale of 1 km 3 water equivalent to test the neutrino 
flux predictions above the TeV that arise in a number of 
models attempting to explain the origin of highest energy 
observed cosmic rays and gamma rays 0j . EeV fluxes are 
difficult to avoid both in the production of the highest en- 
ergy cosmic rays and in the propagation through the cos- 
mic microwave background. Moreover there are allowed 
regions in parameter space for neutrino oscillations || 
which may be best probed with high energy neutrinos 
from cosmological or galactic distances . It is thus de- 
sirable to explore possibilities for alternative neutrino de- 
tection such as horizontal air showers Q or radio pulses 
from high energy showers |^^|. These techniques may 
be advantageous at sufficiently high energies || and can 
provide in any case complementary information relevant 
for flavor identification. 

The detection of coherent radio waves from high en- 
ergy showers has been known since the 60's as an inter- 
esting alternative for detecting ultra high energy showers 
M. These showers develop large excess negative charge 
because the vast majority of the shower particles are in 
the low energy regime dominated by electromagnetic in- 
teractions with the electrons in the target (Compton, 
Bhabha, Moller scattering and electron positron annihila- 
tion) . The excess charge becomes about 20% of the total 
number of electrons and positrons (shower size), which 
is proportional to shower energy. This excess charge ra- 



diates coherently. As long as the wavelength is larger 
than the shower dimensions, the electric field amplitude 
E = \E\ scales with shower energy. The technique has 
been proposed for detecting neutrino interactions in ice 
or sand 0. It has potential advantages such as the rel- 
atively low cost of the detectors (antennae), the large 
attenuation length for radio waves and most importantly 
the fact that information on the excess charge distribu- 
tion can, in principle, be reconstructed from the radiation 
pattern because the radiation is coherent. 

When a particle of charge z travels through a medium 
of refractive index n with velocity v = /3c > c/n 
Cherenkov light is emitted at the Cherenkov angle 9c, 
verifying cos 9c = {Pn)~ l , with a power spectrum given 
by the well known Frank- Tamm result B : 
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with dl — c/3dt, the particle track, and a the fine struc- 
ture constant. This is the standard approximation used 
for most Cherenkov applications for wavelengths orders 
of magnitude smaller than the tracks. 

The frequency band over which Cherenkov radiation 
is emitted can extend well beyond the familiar optical 
band if the medium is transparent. As the radiation 
wavelength becomes comparable to the particle tracks 
the emission from all particles is coherent and the ex- 
cess charge distribution in the shower generates a com- 
plex radiation pattern. It is most convenient to work di- 
rectly with the Fourier transform of the radiated electric 
field, E, which can be directly obtained from Maxwell's 
equations in a dielectric medium pd}^13|. In the Fraun- 
hofer limit (observation distance R much greater than the 
tracklength) the contribution of an infinitesimal particle 
track 51 = vSt is given by: 
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where k, k || R, is the wave vector in the direction of 
observation (R) and l± is the tracklength projected onto 
a plane perpendicular to the observing direction. 

This simple expression displays in a transparent form 
three most important characteristics of such signals: The 
proportionality between the electric field amplitude and 
the tracklength, the fact that in the Cherenkov direction 
(u) — k ■ v = 0) there is no phase factor associated to 
the position along the track direction and the fact that 
radiation is polarized in the direction of l±, that is in the 
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apparent direction of the track as seen from an observer 
located at x. 
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FIG. 1. Electric field angular distribution. 

Recent numerical simulations of radio pulses from 
both electromagnetic and hadronic showers ||,|^@ are 
based on this expression. For energies below 10 PeV 
full simulations are possible fTi) ], The characteristic an- 
gular distributions and frequency spectra are shown in 
Figs. |^^. We can understand most of the pulse character- 
istics by studying the particle distributions in a shower as 
the excess charge follows the electron distribution closely. 
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FIG. 2. Electric field frequency spectrum. 

To a good approximation the pulse is the Fourier trans- 
form of the spatial distribution of the excess charge. For 
many purposes it is sufficient to study the Fourier trans- 
form of the one dimensional distribution (in shower axis 
z) as has been extensively checked Q: 



We here introduce the parameter p(9,u) = u/c(l — 
n cos 9) to transparently relate the radio emission spec- 
trum to the Fourier transform of the (excess) charge dis- 
tribution Q(z). This approximation together with hybrid 
techniques combining simulation and parameterization of 
shower development have allowed the characterization of 
pulses from showers of energy up to 100 EeV |l6|]l7f l . 

The angular distribution of the pulse has a main 
"diffraction" peak corresponding top — 0, the Cherenkov 
direction, see Fig. |[ For \p\ l s h <C 1, where l s h is 
a length scale parameter for the shower fll6f , the elec- 
tric field spectrum accurately scales with electron track- 
length, see Fig. ||. In electromagnetic showers the track- 
length is proportional to the energy |16| and for hadronic 
showers it scales with a slowly varying fraction of the en- 
ergy (80 — 92 % for shower energies between 100 TeV and 
100 EeV) @. 

The scaling with electromagnetic energy is broken by 
interference from different parts of the shower when 
\p\ l s h ~ 1. As a result the frequency spectrum stops 
rising linearly with frequency and has a maximum u>m(9) 
which depends strongly on 9 as seen in Fig. |^. Expand- 
ing the condition for p(9) about 9c it simply becomes 
nsm9cl s hA9uJM / c ~ 1 which clearly displays how A9 
is inversely proportional to ujm as shown in the figure. 
This allows independent establishment of the angle be- 
tween the observation and Cherenkov directions which is 
not sufficient to establish the shower direction but it can 
be combined with other measurements to provide useful 
information. This relation however breaks down when 
approaching the Cherenkov direction because the lateral 
distribution plays the destructive role although there is 
no interference from different shower depths (Eq. ||) . 

The " central peak" at 1 GHz concentrates most of the 
power. For given frequency the angular spread of the 
pulse is also inversely proportional to l s i, . This effect 
hardly shows in showers below 10 PeV with a longitudinal 
scale that only depends logarithmically on energy p9| . 
The difference between the longitudinal development of 
the excess charge for electromagnetic and hadronic show- 
ers is not enough to show up in the radiopulse structure 
(both are governed by the radiation length of the mate- 
rial) . Nevertheless the angular width of the pulse reduces 
significantly for the characteristically elongated electro- 
magnetic showers above 100 PeV because of the LPM 
effect pG[ ]. This narrowing of the angular distribution 
allows the identification of elongated showers. 

The LPM effect manifests as a dramatic reduction of 
the pair production and bremsstrahlung cross sections at 
large energies due to large scale correlations in the atomic 
electric fields It only affects the development of 

showers initiated by photons, electrons or positrons above 
a given energy, about 20 PeV in ice j2Q|. Showers initi- 
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ated by EeV hadrons have high multiplicities (50-100) in 
their first interaction, and the pions produced typically 
have energies 1 — 2% that of the initial hadron. More- 
over 7r°'s above 6.7 PeV are more likely to interact in ice 
than to decay and only about 2% of the hadron showers 
above 10 EeV have one photon with more than 10% of the 
hadron energy. Furthermore in a 100 EeV neutrino inter- 
action for example the fraction of energy transferred to 
the hadron debris (25% in average) fragments into about 
17 hadrons (mostly pions) which have about 5% of the 
transferred energy (except for the leading baryon which 
would carry a fraction X — K where K is the inelasticity). 
As a result the photons that are responsible for the elec- 
tromagnetic subshowers (from ir° and other short lived 
particles decay) have energies which are far removed from 
that of the initial neutrino. Very few hadronic showers 
induced by neutrino interactions of 100 EeV would dis- 
play an LPM tail. For the photon to exceed 100 PeV 
with a probability greater than 2%, the initial neutrino 
energy should be above 80 EeV. 
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FIG. 3. Longitudinal development of electromagnetic and 
hadronic showers. 

The elongation has a dramatic effect on the angu- 
lar distribution of the radio pulse. For electromagnetic 
showers the central peak width narrows as E~ x l z above 
20 PeV |l6[| . A 10 EeV electron produces a pulse which is 
about 10 times narrower than that of a hadronic shower 
of the same energy what makes differentiation between 
pulses from electromagnetic and hadronic cascades pos- 
sible in principle, allowing the characterization of elec- 
tron neutrinos (see Fig. ^). For showers initiated by 
hadrons above 10 EeV the pulse shows a characteristic 
angular distribution of interference of two periodicities 
corresponding to the two length scales, one associated to 
the hadronic shower and the second, longer but of much 
less intensity associated to the electromagnetic LPM tail 
p7[ . The radio pulse for an electron neutrino interaction 



has an interference pattern of similar nature. As the av- 
erage fraction of energy transfer to the hadron debris is 
expected to be about < y >— 0.25 pll this interference 
effect is typically enhanced as shown in Fig. [|. The an- 
gular distribution of the pulse retains enough information 
to allow independent extraction of the total electromag- 
netic energy in both showers, that is to determine the 
individual energy transfer of the reaction. 

Typical energy thresholds for detecting electromag- 
netic showers with single antennas have been made in 
HHJl^g and are typically tens of PeV for showers 
produced at 1 km from the antenna assuming nominal 
frequencies of vq = 1 GHz and band widths of 0.1 vq. 
This corresponds of course to the case that the antenna 
lies just in the illuminated region of the central peak. The 
volume of the illuminated region decreases linearly as v$ 
rises and is also significantly different for electromagnetic 
and hadronic showers of energy above 100 PeV. 

Clearly a signal from a single antenna would be of lit- 
tle use for neutrino detection unless information about 
the shower direction and/or the shower energy could be 
obtained from it. If this was not the case it would be im- 
possible to distinguish them from nearby pulses produced 
by low energy showers such as those induced by deeply 
penetrating muons. Information on neutrino interactions 
can only be obtained by placing antennas in an array cov- 
ering a large region whether on the ice surface or under 
it |l(|. The arrival times for pulses, the polarization [ pT| , 
the relative amplitudes of the signals, and the frequencies 
at different positions of the array elements are in princi- 
ple experimentally accessible and would give relevant and 
redundant information. The technique is similar to "con- 
ventional" neutrino detector proposals but can be highly 
enriched with the angular diffraction patterns, the fre- 
quency spectra and the polarization. 

For intermediate energies one looks for events coming 
from "below" where the Earth provides a shield for all 
other types of particles. For extremely high energies (> 
PeV) however the Earth becomes opaque and neutrino 
events have to be searched in the horizontal direction or 
possibly from " above" . Although some high energy show- 
ers can be expected from other processes such as atmo- 
spheric muon bremsstrahlung at PeV energies the back- 
ground of these events is sufficiently suppressed. There 
is redundant information that allows a variety of cross 
checks. For instance timing can be used to establish the 
shower position in a manner very similar to arrays of 
particle detectors detecting air showers, but the spatial 
distribution of the signal in the detector can be also used 
for the same purpose, even signal polarization provides 
an interesting cross check of the shower orientation, what 
will also be particularly useful to filter spurious noise sig- 
nals out. 
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FIG. 4. Interference pattern for v — 100 MHz as obtained 
in charged current neutrino interactions with different energy 
transfers to the hadron (y) as indicated. 



trino interactions from the rest, and in principle how 
the technique can be used to extract the energy frac- 
tion transmitted to the electron. We have avoided the 
discussion of unresolved experimental issues ]0^3), i.e. 
noise, which are likely to determine the final sensitivity 
of the technique, that is the precise energy value above 
which showers become detectable over sufficiently long 
distances. This sensitivity will be also completely de- 
pendent on the experimental setup which will have to 
be optimized accordingly. These crucial issues have to 
be addressed with in situ experiments and there are ef- 
forts in this direction p5j , but are unlikely to change the 
general conclusions obtained here. 

Acknowledgements: We thank F. Halzen for sug- 
gestions after reading the manuscript and G. Parente, 
T. Stanev and I.M. Zheleznykh for helpful discussions. 
This work was supported in part by CICYT (AEN96- 
1773) and by Xunta de Galicia (XUGA-20604A96). J. A. 
thanks the Xunta de Galicia for financial support. 



The antenna/array parameters are crucial for perfor- 
mance. Most importantly operating frequency vq, band- 
width Av and array spacing. These parameters are 
deeply interconnected at detection level and therefore re- 
quire complex optimization once the noise levels are well 
understood. The spacing of the antenna array will deter- 
mine the minimum distance at which the geometry of the 
illuminated region corresponding to the Cherenkov cen- 
tral peak can be reconstructed, what would indicate the 
position of the shower The nominal frequency will 
determine both the width of the diffraction peaks and 
the transmitting properties of the medium so that ar- 
ray spacing should be adjusted to the choice of antenna. 
Lastly the larger the bandwidth the better the signal to 
noise ratio because the noise behaves as V Aza 

Neglecting attenuation the ratio of shower energy to 
distance has to be above a given value for a shower to 
be detected. As v§ approaches the naively optimal value 
for coherence of 1 GHz the attenuation distance drops 
below 1 km and what is more problematic temperature 
effects become important |ll| possibly forcing detections 
to be within scales less than 1 km for such antennas. 
Using lower frequencies the signal to noise ratio drops 
and higher energy thresholds are needed to compensate 
the loss of signal. This may be nevertheless advisable if 
one is ready to concentrate on neutrinos of the highest 
energies allowing detection at distances above 1 km. 

We have discussed the implications of radio pulse cal- 
culations for high energy shower detection stressing how 
different features of the signal can be used for shower 
characterization. We have shown how the LPM effect 
allows the separation of charged current electron neu- 
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